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O I Abstract. We prove that in the infinite speed-of-light limit (i.e., non-relativistic and 

■ subhorizon limits) , the relativistic fully nonlinear cosmological perturbation equations 
Ci • in two gauge conditions, the zero-shear gauge and the uniform-expansion gauge, exactly 

reproduce the Newtonian hydrodynamic perturbation equations in the cosmological 
background; as a consequence, in the same two gauge conditions, the Newtonian 
> : hydrodynamic equations are exactly recovered in the Minkowsky background. 

p 

■ 1. Introduction 

^ " Newton's gravity and Einstein's gravity are two competing and complementary 

theoretical frames where current cosmological research is practically based. Only with 
i the advent of Einstein's gravity the history of modern physical cosmology has begun [1]. 

H ! Equations describing the cosmological world model were first derived based on Einstein's 

gravity by Friedmann in 1922 [2] , and the Newtonian study followed later by Milne and 
McCrea in 1934 [3]. The two results coincided in the zero-pressure limit. Similarly, the 
linear perturbation equations were first derived based on Einstein's gravity by Lifshitz 
in 1946 [1], and the Newtonian study followed later by Bonnor in 1957 [5]. For density 
perturbation, the two results again coincided in the zero-pressure limit |6]. 

The Newtonian limit of Einstein's gravity in cosmology, however, is not a settled 
subject. The problem is rather serious in the level of background world model where 
the well known Newtonian cosmology [3] is a specially guided one by Einstein's gravity, 
see [?]. Even in the post-Newtonian approach in cosmology we have to subtract the 
background equations based on Einstein's gravity to get proper Newtonian equations 
for the perturbation [8]. 

Despite such a trouble in the background level, curiously we often have quite 
successful relativistic/Newtonian correspondence in the perturbation level. Exact 
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Newtonian perturbation equations are recovered in the Newtonian limit of the post- 
Newtonian approximation independently of the gauge condition [8]; exact Newtonian 
equations are recovered in Minkowsky background [9]. In the perturbation theory, 
however, the analysis and results crucially depend on the gauge choice [lOl [TTl [T2] . 

The exact correspondence of all three Newtonian perturbation variables (perturbed 
density, velocity, and gravitational potential) is not available in a single gauge condition 
even to the linear order [HI [12]. Such a correspondence of all variables is available only 
in the subhorizon limit in the zero-shear gauge and the uniform-expansion gauge |12j . 
which is shown to be valid even to the second order in perturbation |13j . 

However, whether the relativistic/Newtonian correspondence in the subhorizon 
limit is valid to the fully nonlinear order was not known in the literature. Here we 
provide a proof of the correspondence in the two fundamental gauge conditions. For 
the correspondence we will show that we need weak gravity and slow-motion limits, 
in addition to the subhorizon limit as well as negligible pressure and internal energy 
density compared with the energy density; all these limits can be summarized as the 
infinite speed-of-light limit. Such a proof is now available mainly due to a recent advent 
of the exact and fully nonlinear cosmological perturbation theory in Einstein's gravity 

Section |2] is a summary of the fully nonlinear cosmological perturbation theory 
presented in [H]. Section [3] presents our main result, proving the Newtonian limit. 
Section H] is a discussion. 

2. Fully nonlinear perturbations 

We consider the scalar- and vector-type perturbations in a flat background with the 
metric convention [151 E] 

ds^ = - (1 + 2a) c^dt^ - 2xicdtdx' + (1 + 2lp) S^dx'dx^, (1) 

where a{t) is the cosmic scale factor, and a, (p and Xi ^i-re functions of spacetime with 
arbitrary amplitudes; index of Xi is based on 5ij as the metric. The spatial part of 
the metric is simple because we already have taken the spatial gauge condition without 
losing any generality to fully nonlinear order [151 E] ■ 

We consider a fluid without anisotropic stress. The energy momentum tensor is 
given as 

Tab = gC^UaUb + p (gab + UaUb) , (2) 

where tildes indicate covariant quantities; Ua is the normalized fluid four- vector; g and 
p are the mass density and pressure, respectively. We decompose the fluid quantities to 
background and perturbation as 

g = g + 6g, p = p + 6p, Ui = avi, (3) 

where the index of Vi is based on 5ij as the metric. In the explicit presence of the internal 
energy gYi, g should be replaced by ^(1 + H/c^) [9]; in the latter expression g is the 
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material density. We can decompose Xi and Vi to scalar- and vector-type perturbations 
to nonlinear order as 



(4) 



with \1/ 



{v)i 



Xi = X,i + 

= We assign dimensions as 

[a] = [9at] = [ua] = [«] = i:^] = [xi = [v'] = 1, 
[x'] = L, [k] = T-\ [Tab] = [gc^] = [p], [Gg]=T-^, 

where k will be introduced below. 

A complete set of fully nonlinear perturbation equations without taking the 
temporal gauge is 
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(13) 

Equations (l6|)- f|T2|) are the definition of k, energy constraint, momentum constraint, 
trace of ADM propagation, tracefree ADM propagation, covariant energy conservation, 
and covariant momentum conservation equations, respectively. These equations were 
derived in equations (21)-(25), (117), (118), and (26), respectively, of [H]; here we set 
the energy density as /I = qc? and recover c. 

To the background order, equations ([7]), ([9]) and ffTTl) give 



SttG 
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As mentioned, contrary to common beliefs, even with p = = A, these equations are 
pure relativistic in their origin; for a clarification, see |7]. 

Equations (l6l)-(IT2l) are presented without taking the temporal gauge (hypersurface 
or slicing) condition. As the temporal gauge condition we can impose any one of the 
following conditions 



comovmg gauge : 
zero — shear gauge : 
uniform — curvature gauge : 
uniform — expansion gauge 
uniform — density gauge : 



= 0, 
A^O, 
<^ = 0, 
K = 0, 
5 = 0, 



(15) 



or combinations of these to each perturbation order; we may call these the fundamental 
gauge conditions. With the imposition of any of these temporal gauge conditions the 
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remaining perturbation variables are free from the remnant (spatial and temporal) 
gauge mode, and have unique gauge-invariant combinations. Thus, we can regard 
each perturbation variable in those gauges as the gauge-invariant one to nonlinear 
order [151 El]- The temporal synchronous gauge condition sets a = 0; however, as 
this condition fails to fix the gauge degree of freedom completely, and leaves remnant 
gauge mode even from the linear order, we do not consider it as the fundamental gauge 
condition in the nonlinear perturbation analysis. 



3. Newtonian correspondence 

3.1. Newtonian equations 

The Newtonian hydrodynamic equations in the cosmological background are 

^+3-^=--V-(^v), (16) 
a a 

V + -V + -V • Vv = -VU - ^Vp, (17) 
a a a ag 

= -AtvG {g - g) , (18) 

where g is the material density. These equations properly reduce to the Newtonian 
hydrodynamic equations in the Minkowsky background where the background order 
quantities become a = 1 and g = with no A. From equations ( IT6ll and ( IT8l) we have 



f/ + -f/ = — — V - H- 19 
a A 

Our aim in this work is to show that equations (!6|)- (fT2|) properly reduce to equations 
f ll6p -f ll9p in the infinite speed-of-light limit in the zero-shear gauge and in the uniform- 
expansion gauge to the fully nonlinear order. 

Equation ( 1T6|) can be decomposed to the background order part (the third one of 
equation (HM with vanishing pressure term) and the perturbed part 

6=~V-[{l + 6)v], (20) 

where S = Sg/g; Newtonian derivation of equations f lTTj) . f lTS]) and f l2U]) can be found in 
sections 7-9 of [T6l. 



3.2. Infinite speed-of-light limit 

The relativistic perturbation variables [a, ip, w*, g and p) in this paragraph are valid 
in both the zero-shear gauge and the uniform-expansion gauge. As the non-relativistic 
limit we consider 

a < 1, <Cl, < 1, p-C^c^ fl/c^ < 1. (21) 

We identify 

a = -if/, ^=^V, t;^ = iv, (22) 
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where v is the perturbed Newtonian velocity; —U and V correspond to the Newtonian 
and the post-Newtonian perturbed gravitational potentials, respectively [9l [8] ; later we 
will show ip = —a, thus V = U. As the subhorizon limit, we take the dimensionless 
quantity 

^ » 1. (23) 

where H = a/a; in the presence of A we consider ~ SttGq. The first two conditions in 
equation fl2Tl) are the general relativistic (weak gravity) limits, the third one is the special 
relativistic (slow-motion) limit, and the last two conditions imply that we ignore the 
pressure and internal energy density compared with the energy density. Both the non- 
relativistic limits in equation ( pT]) and the subhorizon limit in equation (123!) correspond 
to taking c — > oo limit. 

To the linear order, the vector-type part of equation ([8]) gives 

c^^vj/f) = -SnGgvl'^l (24) 
a'' 

This is "the initial value equation for the frame-dragging potential [10]; '^'f^ is 

supported by v^'"^ which is related to the fluid vorticity [lOj . Thus, we notice that in the 
subhorizon limit "^[^^ is suppressed by a factor in equation f l23|) compared with the vl^^ . 
Therefore, we can ignore ^ in the non- relativistic limit, thus 

X^ = X,- (25) 
while keeping vl;^\ 

3.3. Proof in the Zero-shear gauge 

The zero-shear gauge sets x = [TT]- The condition Xi = leads to vanishing shear 
of the normal-frame vector to the fully nonlinear order [TOl [H], thus we term it as the 
zero-shear gauge; in the literature it is often known as the conformal Newtonian gauge, 
the longitudinal gauge, or the Poisson's gauge, with less clear reasons. 

In the subhorizon limit together with the non-relativistic limits, equations ([8]) and 
(ITOl) . respectively, give 

127rGa , , , 

^ = -^^^ • (H ' (26) 

if = -a. (27) 

Using equations fl2^ and fl27|) and using the Newtonian identifications in equation fl22]) 
we can show the following. Equation ( [TT]) gives equation ( !T6|) . Equation ( IT2|) gives 
equation ( !T7|) . Both equations ([7]) and ([9]) give equation (fT8|) . Finally, equation ([6]) 
gives equation (|T9l) . Therefore, we proved the exact Newtonian correspondence in the 
zero-shear gauge. 
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3.4- Proof in the Uniform- expansion gauge 

The uniform-expansion gauge sets k = 0; this was introduced as the uniform-Hubble 
gauge in [in]. The variable k is perturbed part of the trace of extrinsic curvature, and 
equivalently, perturbed part of the expansion scalar of the normal-frame vector with a 
minus sign. 

In the subhorizon limit together with the non-relativistic limits, equations ([8]) and 
( |T0|) . respectively, give 

127rG'a3 

X = -^5a^V • (H , (28) 
!f = -a. (29) 

Using equations ( 128|) and f l29|) and using the Newtonian identifications in equation ( 122|1 
we can show the following. Equation ( ITTj) gives equation ( !T6|) . Equation ( IT2|) gives 
equation ( IT71) . Both equations (I7j) and ([9]) give equation (fT8l) . Finally, equation ([6]) 
gives equation f|T9|) . Therefore, we proved the exact Newtonian correspondence in the 
uniform-expansion gauge. 

4. Discussion 

In this work we have proved that, in the subhorizon limit together with the non- 
relativistic (weak gravity, slow-motion, and negligible pressure and internal energy 
density compared with the energy density) limits, Newtonian nonlinear perturbation 
equations are exactly recovered from the general relativistic ones in the zero-shear 
gauge and in the uniform-expansion gauge. As a consequence, we also have shown 
that the Newtonian hydrodynamic equations are exactly recovered in the Minkowsky 
background. 

Our proof of the Newtonian limit in the two fundamental gauge conditions shows 
that, in order to have the Newtonian hydrodynamic equations we may set Xi = but 
should have ip = —a, thus the metric should be [9] 

ds'^ = - (l - ^2f/^ c^dt'^ + a'^(l + ^2U^ Sydx'dx^ . (30) 

This may look strange because it is well known that the Newtonian equations can be 
recovered from the Newtonian gravitational potential contained in ^oo only. Whereas, 
our derivation apparently demanded the presence of nontrivial gij as well; the U 
term in g^j is the well-known post-Newtonian correction significantly affecting the light 
deflection, see equation (149) of [8]. Here we would like to make a comment on this 
subtle issue. 

In order to have the conservation equations and the Poisson's equation in f fT6|) - 
( ITSjl from equations ( ITTj) . ( lT2l) and ([9]), we only need the U term in ^oo part; equations 
( ITTj) . ( [T2|) and (Q are the (covariant) energy conservation, the (covariant) momentum 
conservation, and the trace of ADM propagation equations, respectively. Whereas, to 
get equations (JTSl) and (fT9l) from equations ([7]) and IQ, respectively, we require ip = —a 
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which is given by equation (nil|) , thus demanding the presence of U term in gij part which 
is ordinarily the first post-Newtonian contribution; equations ([7]), ([6]) and (fTOj) are the 
energy constraint, the definition of k, and the tracefree ADM propagation equations, 
respectively. 

Therefore, the situation can be summarized as the following. The complete 
Newtonian equations indeed follow from parts of Einstein's equations based on ^oo part 
only. Meanwhile, the validity of the rest of Einstein's equations demands the presence 
of nontrivial gij for the self consistency. 

We point out that in our previous proof, based on the growing mode solutions in the 
zero-pressure medium, of the correspondence to the second order we only have assumed 
the subhorizon limit, thus the correspondence is valid for relativistic (strong gravity and 
fast motion) case [13] , whereas in our present general proof valid to fully nonlinear order 
we additionally have assumed the weak gravity as well as slow-motion limits. 

Our proof of the relativistic/Newtonian correspondence to the fully nonlinear order 
in the two gauges is not necessarily a trivial one. For example, the correspondence is not 
available in other fundamental gauge conditions, like the comoving gauge (f = 0), the 
uniform-curvature gauge = 0), uniform- density gauge (5 = 0), and the synchronous 
gauge (a = 0), even to the linear order in the subhorizon limit [13]: notice that each 
one of these other fundamental gauges sets the fundamental fiuid (density or velocity) 
or potential perturbation variables in Newtonian theory equal to zero. Our result that 
the correspondence is possible only in the two particular gauge conditions is consistent 
with the fact that the post-Newtonian approach is consistent only with the zero-shear 
gauge and the uniform-expansion gauge, see [H [18] . 

For the correspondences we have demanded the non-relativistic as well as the sub- 
horizon limits. As we approach the horizon scale and as soon as the (both special and 
general, including pressure) relativistic effects are not negligible, we have pure Einstein's 
gravity correction terms appearing. In order to properly handle the relativistic effects 
we may have three different but complimentary approaches: the relativistic (nonlinear) 
perturbation theory, the cosmological post-Newtonian approach, and the full-blown 
general relativistic numerical relativity. 

The relativistic/Newtonian correspondence in all cosmological scales was studied 
up to the second-order perturbation in [121 US]- The relativistic/post-Newtonian 
correspondence was studied based on the linear perturbation theory and the first post- 
Newtonian approximation in [18]. The results show that, for the correspondence, 
different gauges are suitable for different variables: e.g., the comoving gauge shows exact 
relativistic/Newtonian correspondences for the density and the velocity perturbations 
up to the second order, whereas the zero-shear gauge shows exact relativistic/Newtonian 
correspondence for the gravitational potential perturbation only to the linear order |13j . 
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